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Q^ ; ABSTRACT 

^\ • Theory suggests the existence of neutrino masses, but little more. Facts are 

^-H , coming close to reveal our fantasy: solar and atmospheric neutrino data strongly 

f~| ' indicate the need for neutrino conversions, while LSND provides an intriguing 

^ , hint. The simplest ways to reconcile these data in terms of neutrino oscilla- 

^^ ■ tions invoke a light sterile neutrino in addition to the three active ones. Out of 

lO , the four neutrinos, two are maximally- mixed and lie at the LSND scale, while 

'"^ ■ the others are at the solar mass scale. These schemes can be distinguished at 

^^ , neutral-current-sensitive solar & atmospheric neutrino experiments. I discuss 

^ ' the simplest theoretical scenarios, where the lightness of the sterile neutrino, the 

OO . nearly maximal atmospheric neutrino mixing, and the generation of Atti^q & 

C*^ ' I^m? atm all follow naturally from the assumed lepton-number symmetry and its 

^^ \ breaking. Although the most likely interpretation of the present data is in terms 

^SZ. • of neutrino-mass-induced oscillations, one still has room for alternative explana- 

Qv , tions, such as flavour changing neutrino interactions, with no need for neutrino 

^\ • mass or mixing. Such flavour violating transitions arise in theories with strictly 

massless neutrinos, and may lead to other sizeable flavour non-conservation ef- 

Q.' fects, such as/i— !'e + 7, /i — e conversion in nuclei, unaccompanied by neutrino- 

I ' less double beta decay. 

Oh! 

(D . 
,-J^ ' 1. Introduction 

>: 

k> , Since the pioneer geochemical experiments of Davis and collaborators, under- 

H I ground experiments have by now provided solid evidence for the solar and the at- 

■ " " ' mospheric neutrino problems, two milestones in the search for physics beyond the 

Standard Model (SM) ElfflttQ. Of particular importance has been the recent confir- 
mation by the Super-Kamiokande collaboration a of the atmospheric neutrino zenith- 
angle-dependent deficit, which has marked a turning point in our understanding of 
neutrinos, providing a strong evidence for z/^ conversions. In addition to the neu- 
trino data from underground experiments there is also some indication for neutrino 
oscillations from the LSND experiment 30. 

Neutrino conversions are naturally expected to take place if neutrinos are massive, 
as expected in most extensions of the Standard Model eI. The preferred theoretical 
origin of neutrino mass is lepton number violation, which typically leads also to lepton 
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flavour violating transitions such as neutrino-less double beta decay, so far unobserved. 
However, lepton flavour violating transitions mav arise unaccompanied by neutrino 
masses 3t3 in models with extra heavy leptons tlHl3 and in supergravity li^. Indeed the 
atmospheric neutrino anomaly can be explained in terms of flavour changing neutrino 
interactions, with no need for neutrino mass or mixing M. Whether or not this 
mechanism will resist the test of time it will still remain as one of the ingredients of the 
final solution, at the moment not required by the data. A possible signature of theories 
leading to FC interactions would be the existence of sizeable flavour non-conservation 
effects, such asyU— ^►e + 7, /i — e conversion in nuclei, unaccompanied by neutrino-less 
double beta decay. In contrast to the intimate relationship between the latter and 
the non-zero Majorana mass of neutrinos due to the Black-Box theorem 113 there is no 
fundamental link between lepton flavour violation and neutrino mass. Barring such 
exotic mechanisms reconciling the LSND (and possibly Hot Dark Matter, see below) 
together with the data on solar and atmospheric neutrinos requires three mass scales. 
The simplest way is to invoke the existence of a light sterile neutrino ElOlla. Out of the 
four neutrinos, two of them lie at the solar neutrino scale and the other two maximally- 
mixed neutrinos are at the HDM/LSND scale. The prototype models proposed in E3'E3 
enlarge the SU{2) ® U{1) Higgs sector in such a way that neutrinos acquire mass 
radiatively, without unification nor seesaw. The LSND scale arises at one-loop, while 
the solar and atmospheric scales come in at the two-loop level, thus accounting for 
the hierarchy. The lightness of the sterile neutrino, the nearly maximal atmospheric 
neutrino mixing, and the generation of the solar and atmospheric neutrino scales all 
result naturally from the assumed lepton-number symmetry and its breaking. Either 
Vf. - Vj. conversions explain the solar data with z/^ - z/^ oscillations accounting for the 
atmospheric deficit tJ, or else the roles of Vr and Ug are reversed tJ. These two basic 
schemes have distinct implications at future solar & atmospheric neutrino experiments 
with good sensitivity to neutral current neutrino interactions. Cosmology can also 
place restrictions on these four-neutrino schemes. 

2. Mechanisms for Neutrino Mass 

Why are neutrino masses so small compared to those of the charged fermionsl 
Because of the fact that neutrinos, being the only electrically neutral elementary 
fermions should most likely be Majorana, the most fundamental kind of fermion. 
In this case the suppression of their mass could be associated to the breaking of 
lepton number symmetry at a very large energy scale within a unification approach, 
which can be implemented in many extensions of the SM. Alternatively, neutrino 
masses could arise from garden-variety weak-scale physics characterized by a scale 
(cr) = O {mz) where {a) denotes a SU{2) ® U{1) singlet vacuum expectation value 
which owes its smallness to the symmetry enhancement which would result if (a) and 
m,, —>■ 0. 



One should realize however that, the physics of neutrinos can be rather different 
in various gauge theories of neutrino mass, and that there is hardly any predictive 
power on masses and mixings, which should not come as a surprise, since the problem 
of mass in general is probably one of the deepest mysteries in present-day physics. 

2.1. Unification or Seesaw Neutrino Masses 

The observed violation of parity in the weak interaction may be a reflection of 
the spontaneous breaking of B-L symmetry in the context of left-right symmetric 
extensions such as the ^^7(2)^ ® SU{2)r ® f/(l)Ei, SU{4) ® SU{2) SU{2) S or 
50(10) gauge groups l21l. In this case the masses of the light neutrinos are obtained 
by diagonalizing the following mass matrix in the basis p, f'^ 



Ml D 
D^ Mn 



:i) 



where D is the standard SU(2) ® U{1) breaking Dirac mass term and Mr = Mj is 
the isosinglet Majorana mass that may arise from a 126 vacuum expectation value 
in S'O(IO). The magnitude of the M^uu term E3 is also suppressed by the left-right 
breaking scale. Ml oc 1/Mr li^. In the seesaw approximation, one finds 

M,eff = ML-DM^'D^. (2) 

As a result one is able to explain naturally the relative smallness of neutrino masses 
since m^, oc 1/Mr. Although M^ is expected to be large, its magnitude heavily 
depends on the model and it may have different possible structures in flavour space 
(so-called textures) c3. In general one can not predict the corresponding light neutrino 
masses and mixings. In fact this freedom has been exploited in model building in order 
to account for an almost degenerate seesaw- induced neutrino mass spectrum c3. 

One virtue of the uniflcation approach is that it may allow one to gain a deeper 
insight into the flavour problem. There have been interesting attempts at formulating 
supersymmetric unifled schemes with flavour symmetries and texture zeros in the 
Yukawa couplings. In this context a challenge is to obtain the large lepton mixing 
now indicated by the atmospheric neutrino data. 

2.2. Weak-Scale Neutrino Masses 

Neutrinos may acquire mass from extra particles with masses O (mz) an therefore 
accessible to present experiments. There is a variety of such mechanisms, in which 
neutrinos acquire mass either at the tree level or radiatively. Let us look at some 
examples, starting with the tree level case. 



2.2.1. Tree-level-generated Neutrino Masses 

Consider the following extension of the lepton sector of the SU{2) ® U{1) theory: 
let us add a set of two 2-component isosinglet neutral fermions, denoted u'^i and Si, 
i = e, /i or r in each generation. In this case one can consider the mass matrix (in 
the basis u, u"^, S) c3 

"0 D " 
D^ M (3) 

M^ /i _ 

The Major ana masses for the neutrinos are determined from 

Ml = DM-^fiM^'^D'^ (4) 

In the limit fi ^ the exact lepton number symmetry is recovered and will keep 
neutrinos strictly massless to all orders in perturbation theory, as in the SM. The 
corresponding texture of the mass matrix has been suggested in various theoretical 
models Ej, such as superstring inspired models tJ. In the latter the zeros arise due to 
the lack of Higgs fields to provide the usual Majorana mass terms. The smallness of 
neutrino mass then follows from the smallness oi fi. The scale characterizing M, unlike 
Mn in the seesaw scheme, can be low. As a result, in contrast to the heavy neutral 
leptons of the seesaw scheme, those of the present model can be light enough as to be 
produced at high energy colliders such as LEP o or at a future Linear Collider. The 
smallness of fi is in turn natural, in t'Hooft's sense, as the symmetry increases when 
/x — ^ 0, i.e. total lepton number is restored. This scheme is a good alternative to the 
smallness of neutrino mass, as it bypasses the need for a large mass scale, present in the 
seesaw unification approach. One can show that, since the matrices D and M are not 
simultaneously diagonal, the leptonic charged current exhibits a non-trivial structure 
that cannot be rotated away, even if we set /i = 0. The phenomenological implication 
of this, otherwise innocuous twist on the SM, is that there is neutrino mixing despite 
the fact that light neutrinos are strictly massless. It follows that flavour and CP are 
violated in the leptonic currents, despite the masslessness of neutrinos. The loop- 
induced lepton flavour and CP non-conservation effects, such as /i ^ e -|- 7 eM, or 
CP asymmetries in lepton-fiavour-violating processes such a.s Z —>■ ef or Z —>■ re^ 
are precisely calculable. The resulting rates may be of experimental interest t 
since they are not constrained by the bounds on neutrino mass, only by those on 
universality, which are relatively poor. In short, this is a conceptually simple and 
phenomenologically rich scheme. 

Another remarkable implication of this model is a new type of resonant neutrino 
conversion mechanism EiJ, which was the first resonant mechanism to be proposed 
after the MSW effect o, in an unsuccessful attempt to bypass the need for neutrino 
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Figure 1: Mechanism for One- loop-induced Neutrino Mass. 
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Figure 2: Mechanism for Two-loop-induced Neutrino Mass 

mass in the resolution of the solar neutrino problem. According to the mechanism, 
massless neutrinos and anti-neutrinos may undergo resonant flavour conversion, under 
certain conditions. Though these do not occur in the Sun, they can be realized in the 
chemical environment of supernovae E3. Recently it has been pointed out how they 
may provide an elegant approach for explaining the observed velocity of pulsars a. 

2.2.2. Radiatively Induced Neutrino Masses 

The prototype one-loop scheme is the one proposed by Zee E3. Supersymmetry 
with explicitly broken R-parity also provides alternative one-loop mechanisms to gen- 
erate neutrino mass arising from scalar quark or scalar lepton exchanges, as shown in 

Fig. m- 

An interesting two-loop scheme to induce neutrino masses was suggested by Babu 
E3, based on the diagram shown in Fig. (^. Note that I have used here a slight 
variant of the original model which incorporates the idea of spontaneous e3, rather 
than explicit lepton number violation. 

Finally, note also that one can combine these mechanisms as building blocks in 
order to provide schemes for massive neutrinos. In particular those in which there 
are not only the three active neutrinos but also one or more light sterile neutrinos. 



such as those in ref. I1£H13. In fact this brings in novel Feynman graph topologies. 

2. 3. Super symmetry: R-parity Violation as the Origin of Neutrino Mass 

This is an interesting mechanism of neutrino mass generation which combines see- 
saw and radiative mechanisms Ea. It invokes supersymmetry with broken R-parity, as 
the origin of neutrino mass and mixings. The simplest way to illustrate the idea is to 
use the bilinear breaking of R-parity SEI in a unified minimal supergravity scheme 
with universal soft breaking parameters (MSUGRA). Contrary to a popular miscon- 
ception, the bilinear violation of R-parity implied by the es term in the superpotential 
is physical, and can not be rotated away e^I. It leads also by a minimization condi- 
tion, to a non-zero sneutrino vev, ^3. It is well-known that in such models of broken 
R-parity the tau neutrino Ut acquires a mass, due to the mixing between neutrinos 
and neutralinos c3. It comes from the matrix 



(5) 



where the first two rows are gauginos, the next two Higgsinos, and the last one denotes 
the tau neutrino. The Vu and Vd are the standard vevs, g's are gauge couplings 
and Ml 2 are the gaugino mass parameters. Since the €3 and the f3 are related, the 
simplest (one-generation) version of this model contains only one extra free parameter 
in addition to those of the MSUGRA model. The universal soft supersymmetry- 
breaking parameters at the unification scale mx are evolved via renormalization group 
equations down to the weak scale O (mz)- This induces an effective non-universality 
of the soft terms at the weak scale which in turn implies a non-zero sneutrino vev f 3 
given as 
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{v'.AM^ + fi'v^AB) (6) 



where the primed quantities refer to a basis in which we eliminate the e^ term from 
the superpotential (but reintroduce it, of course, in other sectors of the theory). 

The scalar soft masses and bilinear mass parameters obey AM^ = and AB = 
at mx- However at the weak scale they are calculable from radiative corrections as 

AM ^ -^m^ln 7 

Stt^ mz 

Note that eq. (^ implies that the R-parity-violating effects induced by v'^ are cal- 
culable in terms of the primordial R-parity-violating parameter €3. It is clear that 
the universality of the soft terms plays a crucial role in the calculability of the fg 
and hence of the resulting neutrino mass E3. Thus eq. (M) represents a new kind of 
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Figure 3: Tau neutrino mass from Broken R-parity vs ^ = (€3^^ + nv^y from ref. 



see-saw scheme in which the M^ of eq. (P is the neutrahno mass, while the role of 
the Dirac entry D is played by the ^3, which is induced radiatively as the parameters 
evolve from mx to the weak scale. Thus we have a hybrid see-saw mechanism, with 
naturally suppressed Majorana z/,- mass induced by the mixing between the weak 
eigenstate tau neutrino and the zino. 

In order to estimate the expected z/^ mass let me first determine the tau neutrino 
mass in the most general supersymmetric model with bilinear breaking of R-parity, 
without imposing universahty of the soft SUSY breaking terms. The z/,- mass depends 
quadratically on an effective parameter ^ defined as ,^ = (e3fd + yUf3)^ oc ^3 character- 
izing the violation of R-parity. The expected rriiy^ values are illustrated in Fig. (^. 
The band shown in the figure is obtained through a scan over the parameter space 
requiring that the supersymmetric particles are not too light. This should be com- 
pared with the cosmologically allowed values of the tau neutrino mass J2 f^u ^ 92Qh'^ 
eV (see below). Note that this only holds if neutrinos are stable. In the present 



22U. 



or 



model the z/^ is expected to decay into 3 neutrinos, via the neutral current 
by slepton exchanges. This decay will reduce the relic z/^ abundance to the required 
level, as long as z/^ is heavier than about 200 KeV or so. Since on the other hand 
primordial Big-Bang nucleosynthesis implies that Uj. is lighter than about an MeV or 
so e3 there is a forbidden gap in this model if the majoron is not introduced. In the 
full version of the model the presence of the majoron allows all neutrino masses to be 
viable cosmologically. 

Back to the simplest model with explicit bilinear breaking of R-parity, note that in 
this model the z/^ mass can be very large. A way to obtain a model with a small and 
calculable Ur mass, as indicated by the simplest interpretation of the atmospheric 
neutrino anomaly in terms of z/^ to Ur oscillations is to assume a SUGRA scheme 
with universality of the soft supersymmetry breaking terms at mx- In this case the 
Ut- mass is theoretically predicted in terms of hi, and can be small in this case due 



to a natural cancellation between the two terms in the parameter ^, which follows 
from the assumed universality of the soft terms at mx- One can verify that m,^^ may 
easily lie in the ten electron-volt range. Lower masses require about two orders 
of magnitude in addition to that which is dictated by the RGE evolution, which 
is certainly not unreasonable. Moreover the solution of the atmospheric neutrino 
anomaly may involve some exotic mechanism, such as the FC interactions tl. 

As a last remark I note that z/g and z/^ remain massless in this approximation. 
They get masses either from scalar loop contributions in Fig. (0) or by mixing with 
singlets in models with spontaneous breaking of R-parity eJ. A detailed study is 
now underway of the loop contributions to z/^ and z/g is underway in Valencia. It 
is important to notice that even when m^^ is small, many of the corresponding R- 
parity violating effects can be sizeable. An obvious example is the fact that the 
lightest neutralino decay will typically decay inside the detector, unlike standard 
R-parity-conserving supersymmetry. This leads to a vastly unexplored plethora of 
phenomenological possibilities in supersymmetric physics e3. 

In conclusion one can see that, of the various attractive schemes for giving neu- 
trinos a mass, only the seesaw scheme requires a large mass scale. It gives a grand 
connection between the very light (the neutrinos) and the very heavy (some unknown 
particles). At this stage is is premature to bet on any mechanism and from this point 
of view neutrinos open the door to a potentially rich phenomenology, since the extra 
particles required have masses at scales that could be accessible to present experi- 
ments. In the simplest versions of these models the neutrino mass arises from the 
explicit violation of lepton number. Their phenomenological potential gets even richer 
if one generalizes the models so as to implement a spontaneous violation scheme. This 
brings me to the next section. 

2.4- Majorons at the Weak-scale 

The generation of neutrino masses will be accompanied by the existence of a 
physical Goldstone boson that we generically call majoron in any model where lepton 
number (or B-L) is an ungauged symmetry which is arranged to break spontaneously. 
Except for the left-right symmetric unification approach, in which B-L is a gauge 
symmetry, in all of the above schemes one can implement the spontaneous viola- 
tion of lepton number. One can also introduce it in a seesaw framework, both with 
SU{2) O f/(l) as well as left-right symmetry §. While in the SU{2) O U{1) case it 
is rather simple, in the case of left-right-symmetric models, one needs to implement a 
spontaneously broken global U(l) symmetry similar to lepton number. One interest- 
ing aspect that emerges in the latter case is that it allows also the left-right scale to 
be relatively low e3. Here I do not consider the seesaw-type majorons, for a discussion 
see ref. a. I will mainly concentrate on weak-scale physics. In all models I consider the 
lepton- number breaks at a scale given by a vacuum expectation value {a) ~ m^eak- 



In all of these models, the weak scale arises as the most natural one and, as already 
mentioned, the neutrino masses when (a) -^ i.e. when the lepton-breaking scale 
vanish eS. 

In any phenomenologically acceptable model one must arrange for the majoron 
to be mainly an SU{2) ® ^(1) singlet, ensuring that it does not affect the invisible 
Z decay width, well-measured at LEP. In models where the majoron has L=2 the 
neutrino mass is proportional to an insertion of (cr), as indicated in Fig. @. In 
the supersymmetric model with broken R-parity the majoron is mainly a singlet 
sneutrino, which has lepton number L=l, so that m,^ oc (a) , where (a) = (j^^), with 
u'^ denoting the singlet sneutrino. The presence of the square, just as in the parameter 
^ in Fig. @, reflects the fact that the neutrino gets a Majorana mass which has 
lepton number L=2. The sneutrino gets a vev at the effective supersymmetry breaking 

bCdie iit'susy ''''weak • 

The weak-scale majorons may have other remarkable phenomenological implica- 
tions, such as the possibility of invisibly decaying Higgs bosons c3 which unfortunately 
I have no time to discuss here (see, for instance e3). 

Finally note that if the majoron acquires a KeV mass (natural in weak-scale 
models) from gravitational effects at the Planck scale c2l it may play obey the main 
requirements to play a role in cosmology as dark matter H. In what follows I will just 
focus on two examples of how the underlying physics of weak-scale majoron models 
can affect neutrino cosmology in an important way. 

2.4.1. Heavy Neutrinos and the Universe Mass 

Neutrinos of mass less than O (100 KeV) or so, are cosmologically stable if they 
have only SM interactions. Their contribution to the present density of the universe 
implies 1^ 

^m,^<92Q,h^eV, (8) 

where the sum is over all isodoublet neutrino species with mass less than O (1 MeV). 
The parameter Q^h"^ < 1, where h"^ measures the uncertainty in the present value of 
the Hubble parameter, 0.4 < h < 1, while fii, = Pu/Pc, measures the fraction of the 
critical density pc in neutrinos. For the i^^ and Uj- this bound is much more stringent 
than the laboratory limits. 

In weak-scale majoron models the generation of neutrino mass is accompanied by 
the existence of a physical majoron, which leads to potentially fast majoron-emitting 
decay channels such as 0'E3 

u' -*u + J . (9) 

as well as new annihilations to majorons, 

u' + v' ^ J + J . (10) 



These could eliminate relic neutrinos and therefore allow neutrinos of higher mass, 
as long as the rates are large enough to allow for an adequate red-shift of the heavy 
neutrino decay and/or annihilation products. While the annihilation involves a diag- 
onal majoron-neutrino coupling g, the decays proceed only via the non-diagonal part 
of the coupling, in the physical mass basis. A careful diagonalization of both mass 
matrix and coupling matrix is essential in order to avoid wild over-estimates of the 
heavy neutrino decay rates, such as that in ref. eB. The point is that, once the neu- 
trino mass matrix is diagonalized, there is a danger of simultaneously diagonalizing 
the majoron couplings to neutrinos. That would be analogous to the GIM mecha- 
nism present in the SM for the couplings of the Higgs to fermions. Models that avoid 
this GIM mechanism in the majoron-neutrino couplings have been proposed, e.g. in 
ref. e3. Many of them are weak-scale majoron models Baa'Ea. A general method to 
determine the majoron couplings to neutrinos and hence the neutrino decay rates in 
any majoron model was first given in ref. e3. For an estimate in the model with 
spontaneously broken R-parity e3 see ref. cj. 

One can summarize that since neutrinos can be short-lived their masses can only 
be really constrained by laboratory experiments based on direct search. The cosmo- 
logical and other bounds are important but require additional theoretical elements in 
their interpretation. 

2.4.2. Heavy Neutrinos and Cosmological Nucleosynthesis 

The number of light neutrino species is restricted by cosmological Big Bang Nucle- 
osynthesis (BBN). Due to its large mass, an MeV stable (lifetime longer than ~ 100 
sec) tau neutrino would be equivalent to several SM massless neutrino species and 
would therefore substantially increase the abundance of primordially produced ele- 
ments, such as ^if e and deuterium ESBE3E3. This can be converted into restrictions 
on the Ur mass. If the bound on the effective number of massless neutrino species is 
taken as iVj, < 3.4 — 3.6, one can rule out z/^ masses above 0.5 MeV e3. If we take 
Nj, < 4.5 E3 the rrij,^ limit loosens accordingly, as seen from Fig. (^, and allows a 
Vr of about an MeV or so. 

In the presence of i/t annihilations the BBN rrty^ bound is substantially weakened 
or eliminated H. In Fig. (^) we also give the expected N^ value for different values of 
the coupling g between i/^'s and J's, expressed in units of 10~^. Comparing with the 
SM g = Q case one sees that for a fixed N^°'^, a wide range of tau neutrino masses 
is allowed for large enough values of g. No z/^ masses below the LEP limit can be 
ruled out, as long as g exceeds a few times 10~^. One can also see from the figure 
that N^ can also be lowered below the canonical SM value N^ = 3 due to the effect 
of the heavy Ur annihilations to majorons. These results may be re-expressed in the 
^^ur ~ 9 plane, as shown in figure |^. We note that the required values of girn^^) fit 
well with the theoretical expectations of many weak-scale majoron models. 
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Figure 4: Effective number of massless SM neutrinos equivalent to tlie lieavy z/,- , as 
given in ref. c3. Non-zero g values (in units of 10~^) can lower N^ with respect to the 
SM case g = (dashed line) due to the effect of annihilations. 
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Figure 5: The BBN- allowed regions for each ]Sf^°-^ lie above the respective curve 
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Figure 6: SSM predictions, from ref. e3 

As we have seen Ur annihilations to majorons may weaken or even ehminate the 
BBN constraint on the tau neutrino mass. Similarly, in some weak-scale majoron 
models decays in eq. (^) may lead to short enough Ur lifetimes that they may also 
play an important role in BBN, again with the possibility of substantially weakening 
or eliminating the BBN constraint on the tau neutrino mass 1^. 

3. Neutrino Indications for New Physics 

The most solid indications in favour of new physics in the neutrino sector come 
from underground experiments on solar and atmospheric neutrinos. Here I will discuss 
mainly the situation as based on published ElOS data, which correspond to a 504 
day solar neutrino data sample H and 535 day atmospheric neutrino data sample S, 
respectively. These were the data first presented at the past Neutrino 98 conference 
in Japan. For more recent information see ref. S and E^. 

3.1. Solar Neutrinos 



The data collected by the Kamiokande, and the radiochemical Homestake, Gallex 
and Sage experiments have no Standard Model explanation. The event rates are 
summarized as: 2.56 ±0.23 SNU (chlorine), 72.2 ±5.6 SNU (Gallex and Sage gallium 
experiments sensitive to the pp neutrinos), and (2.44 ± 0.10) x 10^cm~^s~^ (^B flux 
from Super-Kamiokande) 111. In Fig. (§) one can see the predictions of various stan- 
dard solar models in the plane defined by the ^Be and ®B neutrino fluxes, normalized 
to the predictions of the BP98 solar model 1^. Abbreviations such as BP95, identify 
different solar models, as given in ref. Ell. The rectangular error box gives the 3a error 
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Figure 7: Combined ®B plus hep energy spectrum from ref. E3. 



range of the BP98 fluxes. The values of these fluxes indicated by present data on 
neutrino event rates are also shown by the contours in the figure. The best-fit ^Be 
neutrino flux is negative! Possible non-standard astrophysical solutions are strongly 
constrained by helioseismology studies EaS. Within the standard solar model ap- 
proach, the theoretical predictions clearly lie well away from the 3cr contour, strongly 
suggesting the need for new particle physics in order to account for the data ^. 

The most likely possibility is to assume the existence of neutrino conversions, 
such as could be induced by very small neutrino masses. Possibilities include the 
MSW effect e3, vacuum neutrino oscillations OOO, the Resonant t3t3 Spin-Flavour 
Precession mechanism O and, possibly, flavour changing neutrino interactions E3. 

The 504 day Super-Kamiokande data presents no major surprises, except that 
the recoil energy spectrum produced by solar neutrino interactions was measured for 
the first time c3. The increasing role played rate-independent observables such as the 
spectrum, as well as seasonal and day-night asymmetries, marks a turning point in 
solar neutrino research, which will eventually select the mechanism responsible for 
the explanation of the solar neutrino problem. 

A peculiar feature of the measured spectrum is that it shows more events in the 
highest bins. Barringthe possibly of poorly understood energy resolution effects, 
Bahcall and Krastev E3 have noted that if the flux for neutrinos coming from the 
^He + p — >■ "^He + 6^ + 1^^, the so-called hep reaction is over > 20 times larger than 
the best (but uncertain) theoretical estimates, then this could significantly influence 
the electron energy spectrum produced by solar neutrino interactions in the high recoil 
region, with hardly any effect at lower energies. Fig. |^ shows the ratio of the measured 
to the calculated number of events with electron recoil energy E. The crosses are the 
Super-Kamiokande measurements, while the calculated curves are global fits to the 
data. The horizontal line at 0.37 represents the ratio of the total event rate measured 
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Figure 8: Solar neutrino parameters for 2-flavour active MSW neutrino conversions 
with an enhanced hep flux, as given in ref. 113 



by Super-Kamiokande to the predicted event rate M with no oscillations and only 
^B neutrinos. The total flux of hep neutrinos was varied so as to predict the energy 
spectra in the SM and in various neutrino oscillation scenarios. Comparing with the 
one measured at Super-Kamiokande one sees that the spectra with enhanced hep 
neutrinos provide better fits to the data. However, the required amount seems to 
be just too large 113, though the most recent analysis of 708 days Super-Kamiokande 
spectrum data indicates a somewhat lower hep E3. We look forward to see whether 
this will improve in the next round of data. 

The required solar neutrino parameters Am"^ and sin^ 29 are determined through 
a x^ fif of the experimental data. In Fig. (^) we show the allowed two-flavour regions 
obtained in an MSW fit of the solar neutrino data for the case of active neutrino 
conversions. The data include the chlorine, Gallex, Sage and Super-Kamiokande 
total event rates, energy spectrum, and day-night asymmetry I51l expected in the 
MSW scheme due to regeneration effects at the Earth. The analysis uses the BP98 
model but with an arbitrary hep neutrino flux H. One notices from the analysis that 
rate-independent observables, such as the electron recoil energy spectrum and the 
day-night asymmetry (zenith angle distribution), are playing an increasing role in 
ruling out large regions of parameters. Indeed the most recent 708 days data slightly 
disfavours the small mixing angle (SMA) solution LjO. Another example of an 
observable which has been neglected in most analyses of the MSW effect and which 
could be sizeable for the large mixing angle (LMA) region is the seasonal dependence 
in the solar neutrino flux which would result from the regeneration effect at the Earth 
and which has been discussed in ref. Ej. This should play a more significant role in 
future investigations. 

A theoretical issue which has raised some interest recently is the study of the 
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Figure 9: Solar neutrino survival probability in a noisy Sun, from ref. 



possible effect of random fluctuations in the solar matter density E3E3'Ea. The possible 
existence of noise fluctuations at a few percent level is not excluded by present helio- 
seismology studies. In Fig. (P) we show averaged solar neutrino survival probability 
as a function of Ej ISrv? ^ for sin^ 26' = 0.01. This figure was obtained via a numerical 
integration of the MSW evolution equation in the presence of noise, using the density 
profile in the Sun from BP95 in ref. e3, and assuming that the correlation length Lq 
(which corresponds to the scale of the fluctuation) is Lq = O.IA^, where A^ is the 
neutrino oscillation length in matter. An important assumption in the analysis is that 
^/ree "C -^0 ^ Am, where I free ~ 10 cui is the mean free path of the electrons in the 
solar medium. The fluctuations may strongly affect the ^Be neutrino component of 
the solar neutrino spectrum so that the Borexino experiment should provide an ideal 
test, if sufficiently small errors can be achieved. The potential of Borexino in probing 
the level of solar matter density fluctuations provides an additional motivation for 
the experiment E3. 

The most popular alternative solution to the solar neutrino problem is the vacnnm 
oscillation solution E3 which clearly requires large neutrino mixing and the adjustment 
of the oscillation length so as to coincide roughly with the Earth-Sun distance. This 
solutioniits well with some theoretical models E3 and has been recently re-advocated 
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in ref. O. Fig. |TU| shows the regions of just-so oscillation parameters obtained in a 



recent global fit of the data, including both the rates and the recoil energy spectrum. 
Seasonal effects are expected in this scenario and could potentially be used to further 
constrain the parameters, as described in ref. e3, and also to help discriminating it 
from the MSW scenario. 
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Figure 10: Vacuum oscillation parameters, from ref. Ej 



3.2. Atmospheric Neutrinos 



There has been a long-standing discrepancy between the predicted and measured 
v^ /ue ratio 13 of the fluxes of atmospheric neutrinos E3. The anomaly was found both 
in water Cerenkov experiments, such as Kamiokande, Super-Kamiokande and 1MB a, 
as well as in the iron calorimeter Soudan2 experiment. Negative experiments, such 
as Frejus and Nusex have much larger errors. 

Although individual i/^ or z/g fluxes are only known to within 30% accuracy, the 
^'/x l^e ratio is known to 5%. The most important feature of the atmospheric neutrino 
535-day data sample is that it exhibits a zenith-angle-dependent deficit of muon 
neutrinos which is inconsistent with theoretical expectations. For recent analyses see 
ref. EjH. Experimental biases and uncertainties in the prediction of neutrino fluxes 
and cross sections are unable to explain the data. 

Fig. |1^ shows the measured zenith angle distribution of electron-like and muon- 
like sub-GeV and multi-GeV events, as well as the one predicted in the absence 
of oscillation. It also gives the expected distribution in various neutrino oscillation 
schemes. The thick-solid histogram is the theoretically expected distribution in the 
absence of oscillation, while the predictions for the best-fit points of the various 
oscillation channels is indicated as follows: for v^ — > z/^ (solid line), z/^ -^ z/g (dashed 
line) and z/^ -^ v^ (dotted line). The error displayed in the experimental points is 
only statistical. The analysis used the latest improved calculations of the atmospheric 
neutrino fluxes as a function of zenith angle, including the muon polarization effect 
and took into account a variable neutrino production point Bil. 

Clearly the data are not reproduced by the no-oscillation hypothesis. The most 
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Figure 11: Expected zenith angle distributions for SK electron and muon-like sub- 
GeV and multi-GeV events in the SM (no-oscillation) and for the best-fit points of 
the various oscillation channels, from ref. ^. 



popular way to account for this anomaly is in terms of neutrino oscillations. In Fig. 
( |I2| ) I show the allowed parameters obtained in a global fit of the sub-GeV and multi- 
GeV (vertex-contained) atmospheric neutrino data 1^ including the 535 day Super- 
Kamiokande data, as well as all other experiments combined at 90 (thick solid line) 
and 99 % CL (thin solid line) for each oscillation channel considered. The two lower 
panels Fig. (|T2p differ in the sign of the Am^ which was assumed in the analysis of the 
matter effects in the Earth for the v^ — » Vs oscillations. Though z/^ -^ v^ oscillations 
give a slightly better fit than z/^ -^ v^ oscillations, at present the atmospheric neutrino 
data cannot distinguish between these channels. It is well-known that the neutral-to- 
charged current ratios are important observables in neutrino oscillation phenomenol- 
ogy, which are especially sensitive to the existence of singlet neutrinos, light or heavy 
E3. The atmospheric neutrinos produce isolated neutral pions (vr^-events) mainly in 
neutral current interactions. One may therefore study the ratios of vr^-events and 
the events induced mainly by the charged currents, as recently advocated in ref. E3. 
This minimizes uncertainties related to the original atmospheric neutrino fluxes. In 
fact the Super-Kamiokande collaboration has estimated the double ratio of vr^ over 
e-like events in their sample □ and found R = 0.93 ± 0.07 ± 0.19, consistent both with 
i'^ to Ur or z/^ to I's channels, with a slight preference for the former. The situation 
should improve in the future. We also display in Fig. ([1^) the sensitivity of present 
accelerator and reactor experiments, as well as that expected at future long-baseline 
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Figure 12: Allowed atmospheric oscillation parameters for all experiments including 
the 535-day Super-Kamiokande data, combined at 90 (thick solid line) and 99 % CL 
(thin solid line) for all possible oscillation channels, from ref. e3. In each case the best- 
fit point is denoted by a star and always corresponds to maximal mixing, a feature 
which is well-reproduced by the theoretical predictions of the models proposed in 
ref. OO. The sensitivity of the present accelerator and reactor experiments as well 
as the expectations of upcoming long-baseline experiments is also displayed. 



(LBL) experiments. The first point to note is that the Chooz reactor t3 data excludes 
the region indicated for the z/^ — > Uf. channel when all experiments are combined at 
90% CL. 

From the upper-left panel in Fig. (|12|) one sees that the regions of i/^ — > z/,- os- 
cillation parameters obtained from the atmospheric neutrino data analysis cannot be 
fully tested by the LBL experiments, as presently designed. One might expect that, 
due to the upward shift of the Am^ indicated by the fit for the sterile case (due to 
the effects of matter in the Earth) it would be possible to completely cover the corre- 
sponding region of oscillation parameters. Although this is the case for the MINOS 
disappearance test, in general most of the LBL experiments can not completely probe 
the region of oscillation parameters indicated by the z/^ -^ Ug atmospheric neutrino 
analysis, irrespective of the sign of Am^ assumed. For a discussion of the various po- 
tential tests that can be performed at the future LBL experiments in order to unravel 
the presence of oscillations into sterile channels see ref. o. 

However appealing it may be, the neutrino oscillation interpretation of the atmo- 
spheric neutrino anomaly is at the moment by no means unique. Indeed, the anomaly 
can be well accounted for in terms of flavour changing neutrino interactions, with no 
need for neutrino mass or mixing EJ. Investigations involving upward through going 
muons by Superkamiokande 1^3 as well as other experiments will play an important 
role in discriminating between oscillations and alternative mechanisms to explain the 
sub and multi-GeV atmospheric neutrino data ^. 

3.3. LSND, Dark Matter & Pulsars 



3.3.1. LSND 

The Los Alamos Meson Physics Facility has searched for z/^ -^ z/g oscillations using 
u^ from jj'^ decay at rest S. The z/g's are detected via the reaction u^P —^ e^ n., cor- 
related with a 7 from np -^ d'y (2.2 MeV). The results indicate z/^ -^ v^. oscillations, 
with an oscillation probability of (0.3llo ;Jq±0.05)%, leading to the oscillation param- 
eters shown in Fig. (p^. The shaded regions are the favoured likelihood regions given 
in ref. Q. The curves show the 90 % and 99 % likelihood allowed ranges from LSND, 
and the limits from BNL776, KARMENl, Bugey, CCFR, and NOMAD. A search 
for Vfj_ —>■ z/e oscillations has also been conducted by the LSND collaboration. Using 
z/^ from Tr"*" decay in flight, the z/g appearance is detected via the charged-current 
reaction C(z/e ,e~)X. Two independent analyses are consistent with the above sig- 
nature, after taking into account the events expected from the z/g contamination in 
the beam and the beam-off background. If interpreted as an oscillation signal, the 
observed oscillation probability of 2.6 ± 1.0 ± 0.5 x 10^'^, consistent with the Un —>■ 



Ue oscillation evidence described above. Fig. |1J compares the LSND region with the 
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Figure 13: Allowed LSND oscillation parameters versus competing experiments 1^ 

expected sensitivity from MiniBooNE, which was recently approved to run at Fer- 
milab EJ'lI. A possible confirmation of the LSND anomaly would be a discovery of 
far-reaching implications. 



3.3.2. Dark Matter 



Galaxies as well as the large scale structure in the Universe should arise from the 
gravitational collapse of fluctuations in the expanding universe. They are sensitive to 
the nature of the cosmological dark matter. The observations of cosmic background 
temperature anisotropies on large scales performed bythe COBE satellite a combined 
with cluster-cluster correlation data e.g. from IRAS lB can not be reconciled with the 
simplest COBE-normalized critical-matter-density (i.e., ^2^ = 1) cold dark matter 
(CDM) model, since it leads to too much power on small scales. Adding to CDM 
neutrinos with mass of few eV (a scale similar to the one indicated by the LSND 
experiment Q) corresponding to Vl^ ^ .2. results in an improved fit to data on the 
nearby galaxy and cluster distribution t23. The resulting Cold + Hot Dark Matter 
(CHDM) cosmological model is the most successful Qm = 1 model for structure 
formation, preferred by inflation. However, other recent data have begun to indicate 
a lower value for ^2^, thus weakening the cosmological evidence favouring neutrino 
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mass of a few eV in fiat models with cosmological constant ^a = 1 ~ ^m 
sky maps of the cosmic microwave background radiation (CMBR) with high precision 
at the MAP and PLANCK missions should b ring more light into the nature of the 
dark matter and the possible role of neutrinos c£ll. Another possibility is to consider 



unstable dark matter scenarios clH. For example, an MeV range tau neutrino may 
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Figure 14: Expected sensitivity of the proposed MiniBooNE experiment e3 



provide a viable unstable dark matter scenario 1123 if the v-r decays before the matter 
dominance epoch. Its decay products would add energy to the radiation, thereby 
delaying the time at which the matter and radiation contributions to the energy 
density of the universe become equal. Such delay would allow one to reduce the 
density fluctuations on the smaller scales purely within the standard cold dark matter 
scenario. However, upcoming MAP and PLANCK missions may place limits on 
neutrino stability liSfl and rule out such schemes. 

3.3.3. Pulsar Velocities 

One of the most challenging problems in modern astrophysics is to find a consistent 
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explanation for the high velocity of pulsars. Observations 
range from zero up to 900 km/s with a mean value of 450 ± 50 km/s. An attractive 
possibility is that pulsar motion arises from an asymmetric neutrino emission during 
the supernova explosion. In fact, neutrinos carry more than 99% of the new-born 
proto-neutron star's gravitational binding energy so that even a 1% asymmetry in 
the neutrino emission could generate the observed pulsar velocities. This could in 
principle arise from the interplay between the parity non-conservation present in 
weak interactions with the strong magnetic fields which are expected during a SN 
explosion t2a't23. However, it has recently been noted t23 that no asymmetry in 



neutrino emission can be generated in thermal equilibrium, even in the presence of 
parity violation. This suggests that an alternative mechanism is at work. Several 
neutrino conversion mechanisms in matter have been invoked as a possible engine 
for powering pulsar motion. They all rely on the polarization E£2l of the SN medium 
induced by the strong magnetic fields 10^^ Gauss present during a SN explosion. This 
would affect neutrino propagation properties giving rise to an angular dependence of 
the matter-induced neutrino potentials. This would lead in turn to a deformation 
of the "neutrino-sphere" for, say, tau neutrinos and hence to an anisotropic neutrino 
emission. As a consequence, in the presence of non- vanishing v^ mass and mixing the 
resonance sphere for the v^ — ^t conversions is distorted. If the resonance surface lies 
between the Vr and Ve neutrino spheres, such a distortion would induce a temperature 
anisotropy in the flux of the escaping tau-neutrinos produced by the conversions, 
hence a recoil kick of the proto-neutron star. This mechanism was realized in ref. Iii3 
invoking MSW conversions e3 with rriy^ > 100 eV or so, assuming a negligible z/g 
mass. This is necessary in order for the resonance surface to be located between 
the two neutrino-spheres. It should be noted, however, that such requirement is at 
odds with cosmological bounds on neutrinos masses unless the r-neutrino is unstable. 
On the other hand in ref. ell a realization was proposed in the resonant spin-flavour 
precession scheme (RSFP) O. The magnetic field would not only affect the medium 
properties, but would also induce the spin-flavour precession through its coupling to 
the neutrino transition magnetic moment lB. 



Perhaps the simplest suggestion was proposed in ref. a where the required pul- 
sar velocities would arise from anisotropic neutrino emission induced by resonant 
conversions of massless neutrinos (hence no magnetic moment). 



Raffelt and Janka tl3 have argued, however, that the asymmetric neutrino emission 
effect was overestimated, because the variation of the temperature over the deformed 
neutrino-sphere is not an adequate measure for the anisotropy of the neutrino emis- 
sion. This would invalidate all neutrino conversion mechanisms, leaving the pulsar 
velocity problem without any known viable solution. One potential way out would in- 
voke conversions into sterile neutrinos, since the conversions would take place deeper 
in the star. However, it is too early to tell whether or not it works lii^. 

4. Reconciling the Neutrino Puzzles 

Physics beyond the Standard Model is required in order to explain solar and atmo- 
spheric neutrino data. While neutrino oscillations provide an excellent fit, alternative 
mechanisms are still viable. Thus it is still too early to tell for sure whether neu- 
trino masses and angles are really being determined experimentally. Here we assume 
the standard neutrino oscillation interpretation of the data. While it can easily be 
accommodated in theories of neutrino mass, in general the angles involved are not 
predicted, in particular the maximal mixing indicated by the atmospheric data. It 
is suggestive to consider a theory with bi-maximal mixing of neutrinos O if the so- 
lar neutrino data are explained in terms of the just-so solution. This is not easy 
to reconcile in a predictive quark-lepton unification scheme that relates lepton and 
quark mixing angles, since the latter are known to be small. For recent attempts to 
reconcile solar and atmospheric data in unified models see ref. Ej'EU. The story gets 
more complicated if one wishes to account also for the LSND anomaly and for the hot 
dark matter llaO'EJ. As we have seen the atmospheric neutrino data requires Am^^^ 
which is much larger than the scale Arn^Q which is indicated by the solar neutrino 
data. This implies that with just the three known neutrinos there is no room, unless 
some of the experimental data are discarded. 

4.1. Almost Degenerate Neutrinos 

The only possibility to fit solar, atmospheric and HDM scales in a world with just 
the three known neutrinos is if all of them have nearly the same mass ta, of about 
~ 1.5 eV or so in order to provide the right amount of HDM cM (all three active 
neutrinos contribute to HDM). This can be arranged in the unification approach 
discussed in sec. 2 using the Ml term present in general in seesaw models. With this 
in mind one can construct, e.g. unified 5'O(10) seesaw models where all neutrinos 
lie at the above HDM mass scale (~ 1.5 eV), due to a suitable horizontal symmetry, 
while the parameters Am'^Q & Am'^atm appear as symmetry breaking effects. An 



interesting fact is that the ratio Atti^q / Am'^atm appears as rric^ /mt^ o. There is no 
room in this case to accommodate the LSND anomaly. To what extent this solution 
is theoretically natural has been discussed recently in ref. 
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4-2. Four-Neutrino Models 

The simplest way to incorporate the LSND scale is to invoke a fourth neutrino. 
It must be SU{2) ® U{1) singlet ensuring that it does not affect the invisible Z decay 
width, well-measured at LEP. The sterile neutrino Vs must also be light enough in 
order to participate in the oscillations together with the three active neutrinos. The 
theoretical challenges we have are: 

• to understand what keeps the sterile neutrino light, since the SU{2) ® U{1) 
gauge symmetry would allow it to have a large bare mass 

• to account for the maximal neutrino mixing indicated by the atmospheric data, 
and possibly by the solar 

• to account from first principles for the scales Am^^^, /S.m?Q and ^fn^sNO/HDM 

With this in mind we have formulated the simplest maximally symmetric schemes, 
denoted as (er)(/i s) ta and (es)(/ir) 113, respectively. One should realize that a 
given scheme (mainly the structure of the leptonic charged current) may be realized 
in more than one theoretical model. For example, an alternative to the model in tJ 
was suggested in ref. c3. There have been many attempts to derive the above phe- 
nom enolog ical scenarios from different theoretical assumptions, as has been discussed 
here 
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Although many of the phenomenological features arise also in other models, here 
I concentrate the discussion mainly on the theories developed in ref. HO. These are 
characterized by a very symmetric mass spectrum in which there are two ultra-light 
neutrinos at the solar neutrino scale and two maximally mixed almost degenerate eV- 
mass neutrinos (LSND/HDM scale), split by the atmospheric neutrino scale MO. The 



HDM problem requires the heaviest neutrinos at about 2 eV mass liisl. These scales 



are generated radiatively due to the additional Higgs bosons which are postulated, 
as follows: A^lsnd/hdm arises at one-loop, while Am^j„ and Am^Q are two-loop 
effects. Since these models pre-dated the LSND results, they naturally focussed on 
accounting for the HDM problem, rather than LSND. However, in the meantime 
the evidence for hot dark matter has weakened, whereas LSND came into play. In 
contrast to the HDM problem, the LSND anomaly, if confirmed, would be a more 
convincing indication for the existence of a fourth light neutrino species, considering 
that the HDM may be accounted for in a three neutrino degenerate scenario. 

The models in ll§0 are based only on weak-scale physics. They explain the 
lightness of the sterile neutrino, the large lepton mixing required by the atmospheric 
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Figure 15: {eT){fi s) scheme: u^ - Vr conversions explain the solar neutrino data and 
^^l - ^s oscillations account for the atmospheric deficit, ref. t3. 
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Figure 16: (es)(/ir) scheme: z/g - Vg conversions explain the solar neutrino data and 
v^ - Vt oscillations account for the atmospheric deficit, ref. 113. 



neutrino data, as well as the generation of the mass splittings responsible for solar and 
atmospheric neutrino conversions as natural consequences of the underlying lepton- 
number-like symmetry and its breaking. They are minimal in the sense that they 
add a single SU{2) ® U{1) singlet lepton to the SM. Before breaking the symmetry 
the heaviest neutrinos are exactly degenerate, while the other two are still massless 
IIB. After the global U(l) lepton symmetry breaks the heavier neutrinos split and the 
lighter ones get mass. The models differ according to whether the Ug lies at the dark 
matter scale or at the solar neutrino scale. In the {eT){fi s) scheme the Ug lies at the 
LSND/HDM scale, as illustrated in Fig. ([T5|) while in the alternative (es)(/ir) model. 
Us is at the solar neutrino scale as shown in Fig. (0) tZl. In the {er)(ji s) case the 
atmospheric neutrino puzzle is explained by u^ to Ug oscillations, while in (es) {fir) it 
is explained by z/^ to z/^ oscillations. Correspondingly, the deficit of solar neutrinos is 
explained in the first case by Ue to u^- conversions, while in the second the relevant 
channel is Ue to Ug. 

The presence of additional weakly interacting light particles, such as our light 
sterile neutrino Ug , is constrained by BBN since the Ug would enter into equilibrium 
with the active neutrinos in the early Universe (and therefore would contribute to 
jsjrnax-^ via ueutriuo oscillations tHj, unless Am'^ sin^29 < 3 x 10~^ eV'^ Here Am^ 
denotes the mass-square difference of the active and sterile species and 9 is the vacuum 



mixing angle. However, systematic uncertainties in the BBN bounds still caution 
us not to take them too literally. For example, it has been argued that present 
observations of primordial Helium and deuterium abundances may allow up to Ni, = 
4.5 neutrino species if the baryon to photon ratio is small 1^. Adopting this as a limit, 
clearly both models described above are consistent. Should the BBN constraints get 
tighter 1^ e.g. NJJ"^^ < 3.5 they could rule out the {eT){fi s) model, and leave out 
only the competing scheme as a viable alternative. However the possible role of a 
primordial lepton asymmetry might invalidate this conclusion, for recent work on this 
see ref. Bll. 

The two models would be distinguishable both from the analysis of future solar 
as well as atmospheric neutrino data. For example they may be tested in the SNO 
experiment I1H3 once they measure the solar neutrino flux ($^'") in their neutral 
current data and compare it with the corresponding charged current value (<l>^*"). If 
the solar neutrinos convert to active neutrinos, as in the (er)(yU s) model, then one 
expects ^^'-^/^^'^ ~ .5, whereas in the {es){fiT) scheme {i^e conversion to z/g ), the 
above ratio would be nearly ~ 1. Looking at pion production via the neutral current 
reaction z/,- + A^ — > z/,- + 7r° + A^ in atmospheric data might also help in distinguishing 
between these two possibilities l3, since this reaction is absent in the case of sterile 
neutrinos, but would exist in the (es)(/ir) scheme. 

If light sterile neutrinos indeed exist one can show that they might contribute to 
a cosmic hot dark matter component and to an increased radiation content at the 
epoch of matter-radiation equality. These effects leave their imprint in sky maps of 
the cosmic microwave background radiation (CMBR) and may thus be detectable 
with the precision measurements of the upcoming MAP and PLANCK missions as 
noted recently in ref. 
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4.3. MeV Tau Neutrino 

In ref. tH3 a model was presented where an unstable MeV Majorana tau neutrino 
naturally reconciles the cosmological observations of large and small-scale density 
fluctuations with the cold dark matter picture. The model assumes the spontaneous 
violation of a global lepton number symmetry at the weak scale. The breaking of this 
symmetry generates the cosmologically required decay of the v^. with lifetime t^^ ~ 
10^ — 10^ sec, as well as the masses and oscillations of the three light neutrinos z/g , 
z/^ and Vs which may account for the present solar and atmospheric data, though this 
will have to be checked. One can also verify that the BBN constraints can be satisfied. 

5. In conclusion 

The confirmation of an angle-dependent atmospheric neutrino deficit provides, 
together with the solar neutrino data, a strong evidence for neutrino physics beyond 



the Standard Model. Small neutrino masses provide the simplest, but not unique, 
explanation of the data. If the LSND result stands the test of time, this would 
be a puzzling indication for the existence of a light sterile neutrino. The two most 
attractive schemes to reconcile underground observations with LSND invoke either 
Vf^ - Vj. conversions to explain the solar data, with z/^ - z/^ oscillations accounting for 
the atmospheric deficit, or the other way around. These two basic schemes have 
distinct implications at future solar & atmospheric neutrino experiments. SNO and 
Super-Kamiokande have the potential to distinguish them due to their neutral current 
sensitivity. 

Allowing for alternative explanations of the data from underground experiments 
one can still live with massless non-standard neutrinos or even very heavy neutrinos, 
which may naturally arise in many models. Although cosmological bounds are a fun- 
damental tool to restrict neutrino masses, in many theories heavy neutrinos will either 
decay or annihilate very fast, thereby loosening the cosmological bounds. From this 
point of view, neutrinos can have any mass presently allowed by laboratory experi- 
ments, and it is therefore important to search for manifestations of heavy neutrinos 
at the laboratory in an unbiased way. 

Last but not least, though most of the recent excitement comes from underground 
experiments, one should note that models of neutrino mass may lead to a plethora of 
new signatures which may be accessible also at accelerators, thus illustrating the com- 
plementarity between the two approaches in unravelling the properties of neutrinos 
and probing for signals beyond the Standard Model e3. 

I am grateful to Theocharis Kosmas for the kind hospitality at the loannina. I 
thank all my collaborators, without whom a lot of the common work reported here 
would not have materialized. This work was supported by DGICYT grant PB95-1077 
and by the EEC under the TMR contract ERBFMRX-CT96-0090. 

1. Invited talks by Lande, Gavrin, Kirsten, and Suzuki at Neutrino 98, Takayama, 
Japan 

2. NUSEX Collaboration, M. Aglietta et al, Europhys. Lett. 8, 611 (1989); Frejus 
Collaboration, Ch. Berger et al, Phys. Lett. B227, 489 (1989); 1MB Collaboration, 
D. Casper et al, Phys. Rev. Lett. 66, 2561 (1991); R. Becker-Szendy et al, Phys. 
Rev. D46, 3720 (1992); Kamiokande Collaboration, H. S. Hirata et ai., Phys. Lett. 
B205, 416 (1988) and Phys. Lett. B280, 146 (1992); Kamiokande Collaboration, 
Y. Fukuda et al, Phys. Lett. B335, 237 (1994); Soudan Collaboration, W. W. M 
Allison et al, Phys. Lett. B391, 491 (1997). 

3. T. Kajita, Invited talk at Neutrino 98, Takayama, Japan 

4. C. Yanagisawa, talk at the International Workshop on Physics Beyond The Standard 
Model: from Theory to Experiment, Valencia, (World Scientific, 1998, ISBN 981- 
02-3638-7), |http://fiamenco.uv.es//val97.htmj 



5. Y. Suzuki, Invited talk at Neutrino 98, Takayama, Japan 

6. C. Athanassopoulos, [LSND Collaboration], Phys. Rev. Lett. 75 (1995) 2650; 



Phys. Rev. Lett. 77 (1996) 3082 ; C. Athanassopoulos et al, "Evidence for nu(mu) 
-^ nu(e) neutrino oscillations from LSND," Phys. Rev. Lett. 81, 1774 (1998) 

7. W.C. Louis [LSND Collaboration], "LSND neutrino oscillation results and implica- 
tions," Prog. Part. Nucl. Phys. 40, 151 (1998). 

8. For reviews see J. W. F. Valle, Gauge Theories and the Physics of Neutrino Mass, 
Prog. Part. Nucl. Phys. 26 (1991) 91-171 

9. J. Bernabeu, A. Santamaria, J. Vidal, A. Mendez, J. W. F. Valle, Phys. Lett. B187 
(1987) 303; J. G. Korner, A. Pilaftsis, K. Schilcher, Phys. Lett. B300 (1993) 381 

10. M. C. Gonzalez-Garcia, J. W. F. Valle, Mod. Phys. Lett. A7 (1992) 477; erratum 
Mod. Phys. Lett. A9 (1994) 2569; A. Ilakovac, A. Pilaftsis, Nucl. Phys. B437 
(1995) 491; A. Pilaftsis, Mod. Phys. Lett. A9 (1994) 3595 

11. D. Wyler, L. Wolfenstein, Nucl. Phys. B218 (1983) 205 

12. R. Mohapatra, J. W. F. Valle, Phys. Rev. D34 (1986) 1642; J. W. F. Valle, Nucl. 
Phys. B (Proc. Suppl.) Bll (1989) 118-177 

13. L.J. Hall, V.A. Kostelecky and S. Raby, Nucl. Phys. B267, 41 5 (1986). 

14. M.C. Gonzalez-Garcia et al, Phys. Rev. Lett. 82 (1999) 3202 [hep-ph/9809"53l| . 

15. J. Schechter and J.W. Valle, "Neutrinoless Double Beta Decay In SU(2) X U(l) 
Theories," Phys. Rev. D25, 2951 (1982); for reviews see A. Morales, "Review 
on double beta decay experiments and comparison with theory," hep-ph/980954C ; 
B. Kayser, "Neutrino mass physics," Presented at 28th Rencontres de Moriond: 
Electroweak Interactions and Unified Theories, Les Arcs, France, 13-20 Mar 1993. 

16. J. T. Peltoniemi, D. Tommasini, and J W F Valle, Phys. Lett. B298 (1993) 383 

17. J. T. Peltoniemi, and J W F Valle, Nucl. Phys. B406 (1993) 409 

18. D.O. Caldwell and R.N. Mohapatra, Phys. Rev. D48 (1993) 3259 

19. See, e.g. R. N. Mohapatra and G. Senjanovic, Phys. Rev. D23 (1981) 165. 

20. J.C. Rati, A. Salam. Phys. Rev. D8 (1973) 1240 

21. M Gell-Mann, P Ramond, R. Slansky, in Supergravity, ed. P.van Niewenhuizen 
and D. Freedman (North Holland, 1979); T. Yanagida, in KEK lectures, ed. O. 
Sawada and A. Sugamoto (KEK, 1979); R. N. Mohapatra and G. Senjanovic, 
Phys. Rev. Lett. 44 912 (1980). 

22. J. Schechter and J. W. F. Valle, Phys. Rev. D22 (1980) 2227 

23. S. Lola and J.D. Vergados, "Textures for neutrino mass matrices in gauge theories," 
Prog. Part. Nucl. Phys. 40, 71 (1998); G. Altarelh and F. Ferugho, "Neutrino 
mass textures from oscillations with maximal mixing," Phys. Lett. B439, 112 
(1998) |hep-ph/9807353| . 

24. A. loannissyan, J. W. F. Valle, Phys. Lett. B332 (1994) 93-99; B. Bamert, C.P. 
Burgess, Phys. Lett. B329 (1994) 289; D. Caldwell and R. N. Mohapatra, Phys. 
Rev. D50 (1994) 3477; D. G. Lee and R. N. Mohapatra, Phys. Lett. B329 (1994) 
463; A. S. Joshipura, Zeit. fur Physik C64 (1994) 31 

25. M. C. Gonzalez-Garcia, J. W. F. Valle, Phys. Lett. B216 (1989) 360. 

26. M. Dittmar, M. C. Gonzalez-Garcia, A. Santamaria, J. W. F. Valle, Nucl. Phys. 
B332 (1990) 1; M. C. Gonzalez-Garcia, A. Santamaria, J. W. F. Valle, ibid. B342 
(1990) 108; J . Gluza, J. Maalampi, M. Raidal, M. Zralek, Phys. Lett. B407 
(1997) 45; J. Gluza, M. Zralek; Phys. Rev. D55 (1997) 7030 

27. G. C. Branco, M. N. Rebelo, J. W. F. Valle, Phys. Lett. B225 (1989) 385; N. Rius, 
J. W. F. Valle, Phys. Lett. B246 (1990) 249 

28. M. Dittmar, J. W. F. Valle, contribution to the High Luminosity at LEP working 
group, yellow repor t CERN-91/02, p. 98-103 

29. R. Alemany et. al. |hep-ph/9307252| , published in ECFA/93-151, ed. R. Aleksan, A. 
Ah, p. 191-211 



30. Opal collaboration, Phys. Lett. B254 (1991) 293 and Zeit. fur Physik C67 (1995) 
365; L3 collaboration, Phys. Rep. 236 (1993) 1-146; Phys. Lett. B316 (1993) 
427, Delphi collaboration, Phys. Lett. B359 (1995) 411. 

31. J. W. F. Valle, Phys. Lett. B199 (1987) 432 

32. A.Y. Smirnov and S.P. Mikheev, "Neutrino Oscillations In Matter With Varying 
Density," In *Tignes 1986, Proceedings, '86 massive neutrinos* 355-372; L. Wolfen- 
stein, Phys. Rev. D20 (1979) 2634. 

33. H. N unokawa, Y.Z. Qi an, A. Rossi, J. W. F. Valle, Phys. Rev. D54 (1996) 4356- 
4363, |hep-ph/960530ll 

34. D. Grasso, H. Nu nokawa and J. W. F. Valle, Phys. Rev. Lett. 81, 2412 (1998), 
|astro-ph/9803002 . 



35. A. Zee, Phys. Lett. B93 (1980) 389 

36. K. S. Babu, Phys. Lett. B203 (1988) 132 

37. J. T. Peltoniemi, and J. W. F. Valle, Phys. Lett. B304 (1993) 147 

38. M.A. Diaz, J.C. Romao, and J. W. F. Valle, NucL Phys. B524 23-40 (1998), 
|hep-ph/9706315 . 

39. F. Vissani and A. Yu. Smirnov, N ucl. Phys. B460 , 37-56 (1996); R. Hempfling, 
NucL Phys. B478, 3 (1996), and |hep-ph/9702'4T2| ; H.P. Nilles and N. Polonsky, 
NucL Phys. B484, 33 (1997); B. de Carlos, P.L. White, Phys.Rev. D55 4222-4239 
(1997); E. Nardi, Phys. Rev. D55 (1997) 5772; S. Roy and B. Mukhopadhyaya, 
Phys. Rev. D55, 7020 (1997); A.S. Joshipura and M. Nowakowski, Phys. Rev. 
D 51, 2421 (1995); T. Banks, Y. Grossman, E. Nardi, and Y. Nir, Phys. Rev. 
D 52, 5319 (1995); F.M. Borzumati, Y. Grossman, E. Nardi, Y. Nir, Phys. Lett. 
B 384, 123 (1996); A. Faessler, S. Kovalenko and F. Simkovic, "Bilinear R parity 
violation in neu trinoless double beta decay," Phys. Rev. D58, 055004 (1998), 
hep-ph/9712535| ; M. Carena, S. Pokorski and C.E. Wagner, "Photon signatures for 
low-energy supersymmetry breaking and broken R parity," Phys. Lett. B430, 281 
(1998) |hep-ph/9801251|; M.E . Gomez and K. T amvakis, |hep-ph/9801348 . 



40. M.A. Diaz, |hep-ph/9711435|, [hep-ph/9712213| ; J.C. Romao, |hep-ph/971236^ and 
J. W. F. Valle, |hep-ph/9808292| . 

41. G.G. Ross, J. W. F. Valle. Phys. Lett. 151B 375 (1985); John Ellis, G. Gelmini, C. 
Jarlskog, G.G. Ross, J. W. F. Valle, Phys. Lett. 150B 142 (1985); A. Santamaria, 
J. W. F. Valle, Phys. Lett. 195B 423 (1987). 

42. J. Schechter, J. W. F. Valle, Phys. Rev. D25 (1982) 774 

43. See, e.g. A.D. Dolgov, S . Pastor, and J. W. F. Valle, Phys. Lett. B383 (1996) 
193-198, |hep-ph/9602'233| ; S. Hannestad, J. Madsen, Phys. Rev. Lett. 76 (1996) 
2848-2851, Erratum-ibid. 77 (1996) 5148; J.B. Rehm, G. G. Raffelt, A. Weiss 
Astron. & Astrophys. 327 (1997) 443-452, A. D. Dolgov, S. H. Hansen, D. V. 
Semikoz, Nucl. Phys. B524 (1998) 621-638 

44. J. C. Romao, J. W. F. Valle, Nucl. Phys. B38 1 (1 992) 87-108 

45. For reviews see J. W. F. Valle, |hep-ph/9712277| and |iep-ph/9603307| ; see also M. A. 
Diaz, M. A. Garcia-Jareno, D. A. Restrepo, J. W. F. Valle, Nucl. Phys. B527 (1998) 
44-60 and F. de C ampos, O. J. P. Eboh, J. Rosiek, J. W. F. Valle, Phys. Rev.D55 
(1997) 1316-1325, |hep-ph 9601269| . 

46. Y. Chikashige, R. Mohapatra, R. Peccei, Phys. Rev. Lett. 45 (1980) 1926 

47. E.K. Akh medov, A.S. Josh ipura, S. Ranfone and J. W. F. Valle, Nucl. Phys. B441 
(1995) 61 |hep-ph/9501248| . 



48. A. Joshipura and J. W. F. Valle, Nucl. Phys. B397 (1993) 105; J. C. Romao, F. 
de Campos, and J. W. F. Valle, Phys. Lett. B292 (1992) 329 

49. R. Barbieri, J. Ellis & M. K. Gaillard, Phys. Lett. 90 B, 249 (1980); E. Akhmedov, 



Z. Berezhiani & G. Senjanovic, Phys. Rev. Lett. 69, 3013 (1992). 



50. V. Berezinskii, J. W. F. Valle Phys.Lett B318 360-366,1993, [lhep-ph/9309214 | 

51. E. Kolb, M. Turner, The Early Universe, Addison- Wesley, 1990, and references 
therein 

52. J. W. F. Valle, Phys. Lett. B131 (1983) 87; G. Gelmini, J. W. F. Valle, Phys. 
Lett. B142 (1984) 181; J. W. F. Valle, Phys. Lett. B159 (1985) 49; A. Joshipura, 
S. Rindani, Phys. Rev. D46 (1992) 3000 

53. A. Masiero, J. W. F. Valle, Phys. Lett. B251 (1990) 273; J. C. Romao, C. A. Santos, 
and J. W. F. Valle, Phys. Lett. B288 (199 2) 311; J.C. Rom ao, A. loannisian and 
J. W. F. Valle, Phys. Rev. D55, 427 (1997) |hep-ph/960740T . 



54. R.F. Carswell, MNRAS 268 (1994) LI; A. Songalia, L.L. Cowie, C. Hogan and M. 
Rugers, Nature 368 (1994) 599; D. Tytler and X.M. Fan, Bull. Am. Astr. Soc. 26 
(1994) 1424; D. Tytler, talk at the Texas Symposium, December 1996. 

55. N. Hata et al., Phys. Rev. Lett. 75 (1995) 3977; C.J. Copi, D.N. Schramm and 
M.S. Turner, Science 267 (1995) 192 and Phys. Rev. Lett. 75 (1995) 3981; K. 
A. Olive and G. Steigman, Phys. Lett. B354 (1995) 357-362; 

56. S. Sarkar, Rep. Prog. Phys. 59, 1493 (1996); P. J. Kernan and S. Sarkar, Phys. 
Rev. D 54 (1996) R3681 

57. G. Fiorentini, E. Lisi, S. Sarkar and F.L. Villante, "Quantifying uncertainties in 
primordial nucleosynthesis without Monte Carlo simulations," Phys. Rev. D58, 
063506 (1998); E. Lisi, S. Sarkar and F.L. Villante, "The big bang nucleosynthesis 
limit on iV„e«trmo," Phys. Rev. D59, 123520 (1999) 

58. A.D. Dolgov, S. Pastor , J.C. RomSo and J. W. F. Valle, Nucl. Phys. B496 (1997) 
24-40, |hep-ph/961050^ . 

59. S. Hannestad, Phys. Rev. D57 (1998) 2213-2218; M. Kawasaki, K. Kohri, K. Sato, 
Phys. Lett. B430 (1998) 132-139; for a review see G. Steigman; in Cosmological 
Dark Matter, p. 55, World Scientific, 1994, ISBN 981-02-1879-6; A.D. Dolgov, 
S.H. Hansen, S. Pastor and D.V. Semikoz, Nucl. Phys. B548 (1999) 385-407 

60. Y. Fukuda et al. [Super-Kamiokande Collaboration], "Measurement of the solar 
neutrino energy spectrum using neutrino electron scattering," Phys. Rev. Lett. 82, 



2430 (1999) |hep-ex/9812011 



61. Y. Fukuda et al. [Super-Kamiokande Collaboration], "Constraints on neutrino os- 
cillation parameters from the measurement of day night solar neutrino fluxes at 
SuperKamiokande," Phys. Rev. Lett. 82, 1810 (1999) |hep-ex/9812009 . 

62. Y. Fukuda et al. [Super-Kamiokande Collaboration], "Evidence for oscillation of 
atmospheric neutrinos," Phys. Rev. Lett. 81, 1562 (1998) [hep-ex/9807003| ; see also 
|hep-ex/9803'00^ and |hep-ex/980500^ 

63. M.B. Smy, "Solar neutrinos with SuperKamiokande," hep-ex/9903034 . 

64. K. Inoue, talk at Neutrino Telescopes Workshop, Venice, Feb. 1999. 

65. J. N. Bahcall, |astro-ph/980816^ 

66. J. N. Bahcall, S. Basu and M. H. Pinsonneault, Phys. Lett. B 433 (1998) 1. 

67. (GONG) J. Christensen-Dalsgaard et al., GONG Collaboration, Science 272 (1996) 
1286; (BP95) J. N. Bahcall and M. H. Pinsonneault, Rev. Mod. Phys. 67 (1995) 
781; (KS94) A. Kovetz and G. Shaviv, Astrophys. J. 426 (1994) 787; (CDF94) V. 
Castellani, S. Degl'Innocenti, G. Fiorentini, L.M. Lissia and B. Ricci, Phys. Lett. 
B 324 (1994) 425; (JCD94) J. Christensen-Dalsgaard, Europhys. News 25 (1994) 
71; (SSD94) X. Shi, D.N. Schramm and D.S.P. Dearborn, Phys. Rev. D 50 (1994) 
2414; (DS96) A. Dar and G. Shaviv, Astrophys. J. 468 (1996) 933; (CDF93) V. 
Castellani, S. Degl'Innocenti and G. Fiorentini, Astron. Astrophys. 271 (1993) 601; 
(TCL93) S. Turck-Chieze and I. Lopes, Astrophys. J. 408 (1993) 347; (BPML93) G. 



Berthomieu, J. Provost, P. Morel and Y. Lebreton, Astron. Astrophys. 268 (1993) 
775; (BP92) J.N. Bahcall and M.H. Pinsonneault, Rev. Mod. Phys. 64 (1992) 
885; (SBF90) I.-J. Sackman, A.I. Boothroyd and W.A. Fowler, Astrophys. J. 360 
(1990) 727; (BU88) J.N. Bahcall and R.K. Ulrich, Rev. Mod. Phys. 60 (1988) 297; 
(RVCD96) O. Richard, S. Vauclair, C. Charbonnel and W.A. Dziembowski, Astron. 
Astrophys. 312 (1996) 1000; (CDR97) F. Ciacio, S. Degl'Innocenti and B. Ricci, 
Astron. Astrophys. Suppl. Ser. 123 (1997) 449. 

68. J.N. Bahcall, M.H. Pinsonneault, S. Basu and J. Christensen-Dalsgaard, 
Phys. Rev. Lett. 78 (1997) 171. 

69. J. N. Bahcah, Phys. Lett. B338 (1994) 276; V. Castehani, et al Phys. Lett. B324 
(1994) 245; N. Hata, S. Bludman, and P. Langacker, Phys. Rev. D49 (1994) 3622; 
V. Berezinsky, Comm. on Nucl. and Part. Phys. 21 (1994) 249 

70. S.L. Glashow and L.M.Krauss, "'Just So' Neutrino Oscillations," Phys. Lett. 190B, 
199 (1987); S.L. Glashow, P.J. Kernan and L.M. Krauss, '"Just so' neutrino oscilla- 
tions are back," Phys. Lett. B445, 412 (1999) 

71. V. Barger, S. Pakvasa, T.J. Weiler and K. Whisnant, "Bimaximal mixing of three 
neutrinos," Phys. Lett. B437, 107 (1998), |hep-ph/980638^ ; S. Davidson and 
S.F. King, "Bimaximal neutrino mixing in the MSSM with a single righthanded 
neutrino," Phys. Lett. B445, 191 (1998). 

72. R.N. Mohapatra and J. W. F. Valle, "Solar Neutrino Oscillations From Super- 
strings," Phys. Lett. 177B, 47 (1986). 

73. E.Kh. Akhmedov, Phys. Lett. B213 (1988) 64-68; C. S. Lim and W. Marciano, 
Phys. Rev. D37 (1988) 1368 

74. E.Kh. Akhmedov, The neutrino magnetic moment and time variations of the solar 
neutrino flux, hep-ph/9705451 

75. J. Schechter, J. W. F. Valle, Phys. Rev. D24 1883, (1981), Err. 
ibid.D25 283, (1982). 

76. P.I. Krastev and J.N. Bahcall, "FCNC solutions to the solar neutrino problem," 
|hep-ph/9703267| . 

77. J.N. Bahcall and P.I. Krastev, "Do Hep Neutrinos Affect The Solar Neutrino Energy 
Spectrum?," Phys. Lett. B436, 243 (1998); R. Escribano, J.M. Frere, A. Gevaert 
and D. Monderen, "Boron abundance and solar neutrino spectrum distortion," Phys. 
Lett. B444, 397 (1998). 

78. G. Fiorentini, V. Berezinsky, S. Degl'Innocenti and B. Ricci, "Bounds on hep neu- 
trinos," Phys. Lett. B444, 387 (1998) |astro-ph/9810083i 

79. M.C. Gonzalez-Garcia, P.C. de Holanda, C. Pena-Garay, and J. W. F. Valle, in 
preparation. 

80. J.N. Bahcall, P.I. Krastev and A.Y. Smirnov, "Wher e do we stand w ith solar neu- 
trino oscillations?," Phys. Rev. D58, 096016 (1998), |hep-ph/98072Te . 

81. J.N. Bahcall, P.I. Krastev and A.Y. Smirnov, "Is large mixing angle MSW the 



solution of the solar neutrino problems?," hep-ph/9905220 



82. P.C. de Holanda, C. Pena-Garay, M.C. Gonz alez-Garcia and J. W. F. Valle, "Sea- 
sonal dependence in the solar neutrino flux," hep-ph/9903473| . 

83. A.B. Balantekin, J.M. Fetter and F.N. Loreti, Phys. Rev. D54 (1996) 3941-3951; 
F. N. Loreti and A. B. Balantekin, Phys. Rev. D50 (1994) 4762; F. N. Loreti et 
al, Phys. Rev. D52 (1996) 6664. 

84. H. Nunokawa, A. Rossi, V. Semikoz, J. W. F. Valle Nucl. Phys. B472 (1996) 
495-517 [see also talk at Neutrino 96, |hep-ph/9610526 | 

85. P. Bamert, C.P. Burgess and D. Michaud, "Neutrino propagation through helio- 
seismic waves," Nucl. Phys. B513, 319 (1998); C.P. Burgess, |hep-ph/97li425|; 



86. 

87. 



89. 



90. 



91. 
92. 

93. 

94. 



C.P. Burgess and D. Michaud, "Neutrino propagation in a fluctuating sun," Annals 

Phys. 256, 1 (1997) and |hep-ph/9611368| . 

C. Arpesella at al.. Proposal of the Borexino experiment (1991). 

S.P. Mikheyev, A.Yu. Smir nov Phys.Lett. B 429 (1998) 343-348; B. Paid, G. L. 

Fogli, E. Lisi, D. Montanino, |hep-ph/9805293| 

T. K. Gaisser, F. Halzen and T. Stanev, Phys. Rep. 258, 174 (1995) 

M. C. Gonzalez-Garcia, H. Nunok awa, O. L. G. P eres, T. Stanev, J. W. F. Valle, 

Phys. Rev. D58 (1998) 033004 [jhep-ph 9801368| ; for the 535 days data sample 

update, and the comparison of active and sterile channels see M.G. Gonzalez-Garcia, 

H. Nunokawa, O.L. Peres and J. W. F. Valle, "Active-active and active - sterile 

neutrino oscillation solutions to the atmospheric neutrino anomaly," Nucl. Phys. 

B543, 3 (1998), |hep-ph/9807305i 

R. Foot, R. R. Volkas, O. Yasuda, TMUP-HE L-9801; O. Yasud a, |hep-ph/9804400| ; 

G. L. Fogh, E. Lisi, A. Marrone, G. Sc ioscia, |hep-ph/9808205| ; E.Kh. Akhmedov, 

A.Dighe, P. Lipari and A.Yu. Smirnov, hep-ph/9808270 . 

V. Agrawal et al, Phys. Rev. D53, 1314 (1996); L. V. Volkova, Sov. J. Nucl. Phys. 

31, 784 (1980); T. K. Gaisser and T. Stanev, Phys. Rev. D57 1977 (1998). 

A. Smirnov a nd F. Vissani, Ph ys.Lett. B432 (1998) 376; J G. Learned, S. Pakvasa 

and J. Stone, |hep-ph/9805343| ; H. Murayama and L. Hall, |hep-ph/980621g 



hep-ex/9711002 



95. 



96. 
97. 



99. 



;;i992) L1-L5; E.L. Wright et al., Astrophys. 

1994), ed. A. 



M. Appollonio et al. Phys. Lett. B420 397(1998), 

M. C. Gonzalez-Garcia, talk at International Workshop on Particles in Astrophysics 

and Cosmology: From Theory to Observation, Valencia, Spain, May 3-8, 1999, To 

be published in Nucl. Phys. B (Proc. Suppl.), Ed. V. Berezinsky, G. Raffelt and J. 

W. F. Valle. 

Y. Fukuda et al.. Measurement of the Flux and Zenith Angle Distribution of Upward 

Through Going Muons by Superkamiokande. Phys. Rev. Lett. 82, 2644 (1998), 

hep-ex/9812014 

N. Fornengo, M. C. Gonzalez-Garcia, J. W. F. Valle, in preparation 

W. Louis, private communication 

G. F. Smoot et al., Astrophys. J. 396 

J. 396 (1992) L13 

R. Rowan-Robinson, in Cosmological Dark Matter, (World Scientific, 

Perez, and J. W. F. Valle, p. 7-18, ISBN 981-02-1879-6 

100. J.R. Primack and M.A. Gross, "Cold -|- hot dark matter after SuperKamiokande," 
|astro-ph/98r020^ ; E. Gawiser and J. Silk, |astro-ph/9806T97[ Science, 280, 1405 
(1998), and references therein. 

101. S. Hannestad and G. Raffelt, "Imprint of sterile neutrinos in the cosmic microwave 
background radiation," Phys. Rev. D59, 043001 (1999) |astro-ph/9805223 , 

102. G. Gelmini, D.N. Schramm and J. W. F. Valle, "Majorons: A Simultaneous So- 
lution To The Large And Small Scale Dark Matter Problems," Phys. Lett. 146B, 
311 (1984). 

103. J. Bond and G. Efstathiou, Phys. Lett. B265 (1991) 245; M. Davis et al.. Nature 
356 (1992) 489; S. Dodelson, G. Gyuk and M. Turner, Phys. Rev. Lett. 72 (1994) 
3754; H. Kikuchi and E. Ma, Phys. Rev. D51 (1995) 296; H. B. Kim and J. E. 
Kim, Nucl. Phys. B433 (1995) 421; M. White, G. Gelmini and J. Silk, Phys. Rev. 
D51 (1995) 2669; A. Masiero, D. Montanino and M. Peloso, "Can unstable rehcs 
save pure cold dark matter?," [hep-ph/9902380| . 

104. S. Hannestad, "Probing neutrino decays with the cosmic microwave background," 
Phys. Rev. D59, 125020 (1999) |astro-ph/9903475i 

105. A.G. Lyne and D.R. Lorimer, Nature 369 (1994) 127. 



106. N.N. Chugai, Pis'ma Astron. Zh.lO, 87 (1984). 

107. A. Vilenki n, Astrophys. J. 451 (1995) 700; Dong Lai, Y.-Z. Qian, |astro- 
ph/971204^ 



108. A. Kusenko, G. Segre and A. Vilenkin, "Neutrino transport: No asymmetry in 
equilibrium," Phys. Lett. B437, 359 (1998) |astro-ph/9806205i 

109. H. Nunokawa, V.B. Semikoz, A.Yu. Smirnov and J. W. F. Valle, Nucl. Phys. 
B501 (1997) 17 

110. A. Kusenko, G. Segre, Phys. Rev. Lett. 77 (1996) 4872 & Phys. Rev. Lett. 79 
(1997) 2751; Y.Z. Qian, Phys. Rev. Lett. 79 (1997) 2750 

111. E.Kh. Akhmedov, A. Lanza and D.W. Sciama, Phys. Rev. D56 (1997) 6117 

112. H.T. Janka and G.G. Raffelt, "No pulsar kicks from deformed neutrinospheres," 
Phys. Rev. D59, 023005 (1999) |astro-ph/9808099l 



113. D. Grasso, H. Nunokawa, A . Rossi, A. Yu. Smirnov, J. W. F. Valle, in preparation 

114. G. Altarelli and F. Feruglio, Phys. Lett. B451 (1999) 388 |hep-ph/9812475| ; S. Lola 
and G.G. Ross, "Neutrino masses from U(l) symmetries and the SuperKamiokande 
data," hep-ph/9902283 ; R. Barbieri, L.J. Hall and A. Strumia, "Textures for at- 
mospheric and solar neutrino oscillations," Phys. Lett. B445, 407 (1999), |hep 



ph/9808333| ; M.E. Gomez, G.K. Leontaris, S. Lola and J.D. Ver gados, "U(l) tex- 
tures and lepton flavor violation," Phys. Rev. D59, 116009 (1999), p^ph/9810291| ; 
G.K. Leontaris, S. Lola, C. Scheich and J.D. Vergados, "Textures for neutrino mass 
matrices," Phys. Rev. D53, 6381 (1996) 

115. J. A. Casas, J.R. Espinosa, A. Ibarra and I. Navarro, hep-ph/9905381 ; J. Ellis and 
S. Lola, "Can neutrinos be degenerate in mass?," |hep-ph/990427S . 

116. Q. Y. Liu, A. Yu. Smirnov, Nucl. Phys. B524 (1998) 505-523; V. Barger, K. 
Whisnant and T. Weiler, Phys.Lett. B427 (1998) 97-104; S. Gibbons, R. N. Moha- 
patra, S. Nandi and A. Raichoudhuri, Phys. Lett. B430 (1998) 296-302; Nucl.Phys. 
B524 (1998) 505-523; S. Bilenky, C. Giunti and W. Grimus, Eur. Phys. J. C 1, 
247 (1998); S. Goswami, Phys. Rev. D 55, 2931 (1997); N. Okada and O. Yasuda, 
Lit. J. Mod. Phys. A12 (1997) 3669-3694 

117. E. J. Chun, A. Joshipura and A. Smirnov, in Elementary Particle Physics: Present 
and Future (World Scientific, 1996), ISBN 981-02-2554-7; P. Langacker, "A Mech- 
anism for ordinary sterile neutrino mixing," Phys. Rev. D58, 093017 (1998); M. 
Bando and K. Yoshioka, Prog. Theor. Phys. 100, 1239 (1998) 

118. J. R. Primack, et al. Phys. Rev. Lett. 74 (1995) 2160 

119. J. Schechter and J. W. F. Valle, Phys. Rev. D21 (1980) 309 

120. R. Barbieri and A. Dolgov, Phys. Lett. B 237, 440 (1990); K. Enquist, K. 
Kainulaine n and J. Maalam pi, Phys. Lett. B 249, 531 (1992); D. P. Kirilova and 
M. Chizov, |hep-ph/970728^ . 

121. R. Foot, R.R. Volkas, Phys. Rev. D55 (1997) 5147-5176 

122. SNO collaboration, E. Norman et al. Proc. of The Fermilab Conference: DPF 92 
ed. C. Albright, P. H. Kasper, R. Raja and J. Yoh (World Scientific), p. 1450. 

123. A. Joshipura, J. W. F. Valle, Nucl. Phys. B440 (1995) 647. 



